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wrv, Abstract. We investigate the effect of an important class of glitches occurring in 

I— I' the detector data on matched filter searches of gravitational waves from coalescing 

, ■ compact binaries in the advanced detector era. The glitches, which can be modeled 

^ , as sine-Gaussians, can produce triggers with significant time delays and thus have 

OO ' important bearing on veto procedures as will be described in the paper. We provide 

3^ , approximated analytical estimates of the trigger SNR and time as a function of the 

f— ^ I parameters describing the sine-Gaussian (center time, center frequency and Q-factor) 

l' ' and the inspiral waveform (chirp mass). We validate our analytical predictions through 

f~^ , simple numerical simulations, performed by filtering noiseless sine-Gaussians with the 

CO ' inspiral matched filter and recovering the time and value of the maximum of the 

". . ■ resulting SNR time series. Although we identify regions of the parameter space in which 

^ I each approximation no longer reproduces the numerical results, the approximations 

complement each other and together effectively cover the whole parameter space. 
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1. Introduction 

Advanced interferometric gravitational-wave detectors are currently under construction 
and are expected to start delivering strain data in a few years. Major improvements over 
the first-generation interferometers will include a significantly better sensitivity curve at 
all frequencies of interest, which will effectively shift the sensitive band lower boundary 
from ~ 40 Hz to ~ 10 Hz [H [2]. 

A promising source of gravitational- wave detections in the advanced interferometer 
era is represented by the inspiral and coalescence of compact binary systems (neutron 
star binaries, black hole-neutron star binaries and black hole binaries) [3]. The strain 
waveform associated with such events can be calculated to a good approximation at least 
for some regions of the parameter space describing the binary system |1]. Consequently, 
the classical data analysis pipeline for detecting such signals employs a filter matched 
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to the analytical waveform [3]. This approach consists in correlating the strain data 
with the waveform and producing an estimate of the signal to noise ratio (SNR) for 
that waveform. Triggers are then generated whenever the SNR exceeds a pre-defined 
threshold. Lists of triggers are independently produced for each detector and combined 
to obtain coincident triggers, which are then checked for consistency and passed on 
to more elaborate analysis. Similar search pipelines will be employed in the advanced 
detector era. However, in order to exploit the improved noise curves and extend the 
search to lower frequencies, the strain data will have to be filtered with significantly 
longer inspiral waveforms than used for past searches (up to tens of minutes as opposed 
to about one minute for BNS). In fact, the duration of the waveform in the interferometer 
output depends strongly on the low- frequency cutoff [6]. 

It is well known that first generation interferometers were affected by high rates 
of short- duration transients unrelated to astrophysical events and commonly known 
as glitches, most of which can be roughly modeled as short sine-Gaussian waveforms 
with a few cycles, i.e. a small Q-factor [7]. The typical effect of glitches on searches 
for binary coalescence was a significant deviation of the SNR distribution from what 
one would expect for perfectly Gaussian noise, i.e. a x^ distribution. Most glitches 
were ruled out as gravitational-wave candidates by monitoring auxiliary channels and 
their origin could rather be attributed either to the noisy environment of the detector 
(for instance seismic noise, weather, electromagnetic disturbances) or to unexpected 
behavior of the detector components (ADC saturation, problems with control systems 
and thermal compensation, scattered light) [HIE]. Nevertheless, the surviving glitches 
affected the coincidence rate and thus weakened the ability to both confidently detect 
signals as well as set strong upper limits on the coalescence rate. 

Unfortunately, despite the improved sensitivity, advanced detectors are expected 
to manifest similar glitches in their output, which will still pose significant challenges 
to future gravitational-wave searches. In particular, the longer BNS inspiral waveforms 
used to exploit the sensitivity at low frequency could lead to significant time delays 
between the occurrence of a glitch and the coalescence time of the resulting spurious 
trigger. Often, such delays will be much longer than the glitch duration and they will 
have to be taken into account when associating inspiral triggers with excitations of 
the auxiliary channels, or conversely, when using auxiliary channels to veto intervals of 
strain data. 

For this reason, we investigate triggers generated by sine-Gaussian glitches 
interacting with a matched-filter search for inspiral gravitational waves in the advanced 
detector era. To this end we present approximations to the complex matched filter 
output and to the SNR time series, which is the quantity used for producing the inspiral 
triggers. We present three such approximations which complement each other and 
completely cover the parameter space of the sine-Gaussians, namely the central time to, 
the central frequency /o and the quality factor Q. We support the approximations with 
numerical simulations. 

The paper is organized as follows. In section 2 we briefly review the matched 
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filter algorithm and fix the relevant notation. In section 3 we derive the approximated 
estimates for the SNR time series. In section 4 we compare these results to numerical 
simulations and we identify regimes where the different approximations are no longer 
acceptable. Section 5 discusses the implications of our results and possible ways to 
incorporate our estimations into veto procedures for the advanced era. 

2. Matched filter 

A matched filter search combines a detector strain data segment s{t), the signal to be 
sought h(t) and the one-sided power spectrum of the detector noise S^if) into a quantity 
representing the SNR as a function of the parameters describing the gravitational-wave 
source [6]. The SNR for a binary coalescing at time t is calculated aslj 

Pit) = ^ (1) 

where z{t) is the output of the complex matched filter, 

and a is the sensitivity of the detector to the sought after signal, 
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In the above expressions, s{f) is the Fourier transform of the pre-conditioned data 
segment and 

kf) = hof-"'e-'^^f^ (4) 

is the inspiral signal template expressed in the frequency domain via the stationary 
phase approximation [TU]. /l and /h are frequency cutoff values. Searches performed 
on first-generation detector data used /l = 40 Hz and for advanced interferometers 
/l = 10 Hz is expected, /h is usually set to the frequency corresponding to the 
innermost stable circular orbit of a test mass orbiting a Schwarzschild black hole, 
i.e. /isco=cV(6y6vrGM). 

Candidate inspiral triggers are identified as local maxima of the SNR crossing a 
pre-established threshold. Each trigger carries the set of associated parameter values, 
e.g. coalescence time and phase and binary component masses. 

3. Response to a sine-Gaussian waveform 

We are interested in calculating the response of the inspiral matched filter to a sine- 
Gaussian waveform, i.e. the SNR time series p(t) after the center time of the sine- 
Gaussian. From pit) we want to estimate the maximum SNR p and the corresponding 
time t. 

X For simplicity and since search algorithms work with SNR time series, we only explicitly display the 
dependency on the coalescence time t and omit the other parameters. 
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The sine-Gaussian waveform can be written in the time domain as 

s{t) = Aexp (-^t_!^\ cos(27r/ot + 0o) (5) 

where A is an overall amplitude, to and /o are the location of the sine-Gaussian in time 
and frequency respectively, r = Q/27ifQ is the time duration, Q is the dimensionless 
quality factor and 0o is the phase of the oscillation at t = 0. We are free to define the 
time origin as to = 0, so that t represents the time delay between the sine-Gaussian and 
the inspiral trigger. Since the SNR is linear in the signal amplitude, we can also set 
A = 1. The sine-Gaussian can then be expressed in the frequency domain as 



5(/) = ^rexp {-n'r'if - f,f + z0o) 
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exp \-Q ^ - 2z(; 
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For fh < f < fu, the second exponential falls off very quickly as Q increases from 1, so 
we neglect it altogether. We can now also set 00 = since the SNR is not affected by 
global phase factors. Inserting (|4]) and the so-modified ([6]) into (j2]) we get 

z{t) = 2^rh, j ^^ L__exp{-n\\f-fof + zi'{f))df. (7) 

In order to get an explicit expression, we approximate ip{f) with the simple Newtonian 
chirp [H] 

where t is the coalescence time, is the orbital phase at coalescence, ( = ttGA4./c^ 
and A4. = {mim2Y^^ /{nii + 7722)^''^ is the chirp mass of the binary system composed 
of masses mi and m2. Although the resulting integral can not be evaluated exactly, 
different approximations can be found, as will be discussed in the next subsections. 

3.1. Approximation I 

The first approximation we describe makes use of the stationary phase approximation. 
If dip/df = at a frequency /s, we can approximate ipif) around /s with a second order 
power series, 

V^(/)^V'(/s) + /3(/-/s)'. (9) 

For the Newtonian chirp this holds and we have 

P = ys'""C"'' (11) 



where Tq is the first chirp time evaluated using /o as the fiducial frequency [12] • Note 
that /s is a function of time and we are evaluating a family of integrals parametrized by 
t. The factor f~'^'^S^^{f) in (JTj) varies more slowly with frequency than the exponential, 
so we take it as constant and evaluate it at /g. The integrand can then be rewritten as 
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the product of two Gaussians, one real and one complex, so the matched filter output 
is (omitting irrelevant phase factors) 

/.-'/«,/„ (, (j-f,f (f-f„f \ 

^'" " '^^uK) L '"p ['^^f — ^ "^ <'" 



where o"sg = 1/\/2ttt and cTf = y 3/o/167rro(/s//o)'^^^^ are the standard deviations of 
the Gaussians. Using ( TT2l) and (E]) in ([T]), and carrying out the integral using Cauchy's 
theorem, we get 

Pit) ~ TT^e ^^^^ (13) 

withal = J,^f~y^S~\f)df. 

This is a complicated function of t through /g and 0"f. In particular note that t 
also enters via Sn{fs) and so the shape of the SNR pulse depends on the PSD of the 
noise. Although «S'n(/) can in principle be modeled as a combination of power laws, 
the resulting expression is complicated and a numerical maximization is thus needed in 
order to determine p and t. 

3.2. Approximation II 

Our second method is a slightly different application of the stationary phase 
approximation and is expected to be appropriate for small Q. In fact, we still use 
(jni) but we now assume that the largest contribution to the integral in (ITj) comes from a 
narrow band around /g. Within this band, we can thus consider all other terms of the 
integrand constant and equal to their value at /s, including the sine-Gaussian peak, so 
the integral reduces to 

/%W-/.)^d/^,/|. (14) 

Therefore, omitting all irrelevant phase terms, 

,(i) « 2.rfe„ °"P'-f"''^'-^°>'> (15) 

/."°S„(/.)/3V2 

and the resulting SNR time series is 

For the Newtonian chirp, 

nm ~ V6vr3/^rC^/^^exp(-7rV(/.-/o)^) ,, 

This is still a quite complicated function of time and it still depends on the PSD of 
the noise. In general, at least for a smooth PSD with no peaks, it is the product of 
two peaks with different time scales, one corresponding to the exponential term and 
the other associated with the /J''^®5'n(/s)/3^''^ term. Thus i depends on /q, r, ( and 
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the details of S^if) even for the Newtonian chirp. Although we expect approximation 
II to be valid for small Q, we can try to take the result to the limit of Q — )■ oo. The 
exponential becomes more peaked than the f^'^^^S~^{fs)(3~^^'^ term, so p(t) is maximum 
where /s = /o and 

P-^-fo'^'s;;U)r'^' (18) 

i ^ To. (19) 

This suggests that in the large Q limit the trigger delay is just the time the inspiral 
takes to coalesce after crossing the center frequency of the sine-Gaussian, which can be 
expected from intuition. As we said, however, this approximation is expected to work 
for small Q and thus in general p and t must be found numerically. 

3.3. Approximation III 

Although probably less useful in practice, a different approximation can be derived in 
the large Q limit. In this regime the integrand in ([7]) is no longer dominated by the 
band around fg but by the narrow peak of the sine-Gaussian, centered on /q. Again, 
the factor f~'^^^S~^{f) varies slowly with frequency and can be regarded as constant, 
and the integral limits can be extended over the real axis. We still approximate the 
template phase as a second-order Taylor series, but this time around /o, 

^(/) ^ ^ifo) + «(/ - /o) + /?(/ - /o)'. (20) 

For instance, for the Newtonian chirp we have 

a = 2nt- Y^/o"'^\-'/' = Mt - To) (21) 

/3 = ^/o""^'C-^/^ (22) 



After switching the integration variable to x = a/tt^t^ — i(3{f — /o), and omitting 
irrelevant phase factors, we are left with a Gaussian integral 

.W«^#^#^rexpP + ^^)d. (23) 

which readily gives 

27rrho fo''' ( a' \ 

'^^^v/^FV^^^5„(/o)'^Pl, 4(7rV2-./3);- ^'^^ 

The final expression for the approximated SNR time series is 

1 expf-aV4vrV2(l + 4^)) 
p{t) - -— ^-^ ^ ^^- (25) 

-0 fli^s^u.)(^^^^' 



In the case of the Newtonian chirp a is the only function of time, so p{t) is a simple 
Gaussian pulse and we obtain explicit expressions for trigger SNR and time, 

P^^o"Vo"'/'^n^(/o)(l + ^) (26) 

i ^ To. (27) 
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Note that t is consistent with approximation II taken in the limit of large Q, but p is 
not and we expect approximation III to give the correct result in this case. 

4. Numerical simulations 

We test the accuracy of the approximations by numerically computing the response of an 
inspiral matched filter to noiseless sine-Gaussians injected in the time domain. For S'n(/) 
we take the noise PSD corresponding to advanced LIGO's zero detuning, high power 
configuration, an approximation of which is provided by the LALSimulation module of 
LALSuite [13]. For simplicity, and in order to work well within the numerical range of 
the floating-point representation, we set the amplitude of the simulated sine-Gaussian to 
1 and at the same time we scale the noise PSD by 10^^; therefore, the resulting trigger 
SNR refers to sine-Gaussians with an amplitude oi A = 10~^^. We scan the (/o, Q) 
parameter space in the region 10 Hz < /o < min (1 kHz, /isco), 1 < <5 < 1000 for total 
mass values of 2.8Mq, 5Mq and IOMq. For each point, we compute the maximum 
SNR p and the delay t between the maximum SNR and the center of the injected sine- 
Gaussian, and we compare to the predicted values. Results for approximation I, II and 
HI can be seen in figure [H |2] and [3] respectively. 

Approximation I works well for both p and t, but it degrades for large values of /o, 
Q or mass. In particular, we find that surfaces of constant accuracy roughly match those 
of constant f^QM; the 5% accuracy for SNR, at least in the explored parameter range, 
is at f^QM/MQ ^ 5 x 10'^ s~^. This is likely not a major problem, as we expect most 
glitches to last at most tens of cycles and affect mostly low frequencies. Moreover, at 
high frequency or mass the time delay is small and thus the problem we are considering 
is less important. Also note that neglecting the second exponential of (l6l) produces no 
noticeable effects even for Q ~ 1. 

Approximation II produces excellent estimates of the time delay across all the 
explored parameter space, but wrong estimates of SNR for large Q, as can be expected. 
In fact, as Q increases, the sine-Gaussian peak in the frequency domain shrinks and at 
some point the integral is no longer dominated by the region around fg. Since we are 
interested in the low-Q region, despite this problem this is still a useful approximation, 
although not significantly simpler than I. 

As expected, approximation HI works very well in the high mass, high Q region 
where the other two approximations do not give such good results. Moreover it is 
analytically simpler. But for a large region of the parameter space where Q is small, 
this approximation fails and so is not so useful for a large part of the parameter space 
of interest. This is because the width of the Gaussian in the frequency domain becomes 
comparable or larger than the width of the other terms in the integrand and the 
approximation is not valid. However, in this part of the parameter space we can simply 
revert to approximations I and II. 

Although our calculations are based on the Newtonian chirp, search pipelines use 
inspiral waveforms with higher post-Newtonian order [6]. Thus we also check the 
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Figure 1. Comparison between simulations (shaded bands, color online) and 
approximation I (black contours) for different total masses. Left plots show the trigger 
SNR, right ones the delay. The dashed line is the ISCO frequency. As can be seen, 
this approximation fails for large /o, Q or M. 
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Figure 2. Comparison between simulations (shaded bands, color online) and 
approximation II (black contours) for different total masses. Left plots show the trigger 
SNR, right ones the delay. The dashed line is the ISCO frequency. The SNR predicted 
by approximation II is not correct for Q ^ 1, as expected. Remarkably however, the 
predicted trigger time is good for most of the explored parameter space. 
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Figure 3. Comparison between simulations (shaded bands, color online) and 
approximation III (black contours) for different total masses. Left plots show the 
trigger SNR, right ones the delay. The dashed line is the ISCO frequency. As expected, 
approximation III does not work for small Q. 
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Figure 4. Comparison between a simulation using a 3.5PN equal-mass inspiral 
waveform (shaded bands, color online) and approximation I (black contours) for a 
total mass of 5Mq . The left plot is the trigger SNR, the right one is the trigger delay. 
The dashed line is the ISCO frequency. 



accuracy of the approximations against simulations using a 3.5PN inspiral filter. We 
find that the accuracy degrades but is still within a few percent both for p and t and 
it retains a similar dependency on /o, Q and M. An example for a 5Mq equal- mass 
binary is shown in figure HI Since our approximations are based simply on power-series 
approximations of the inspiral phase, they can in principle accommodate high post- 
Newtonian order waveforms, at the price of more complicated expressions for p{t). 

In a search pipeline, the parameter space of the binary is covered by a template 
bank. A true inspiral signal produces triggers only for templates whose parameters 
are close enough to the true values, depending on the ambiguity function of the 
waveform. However, a strong glitch generally excites a significant fraction of the whole 
bank, producing a cluster of triggers with different SNRs and times. An example of 
this phenomenon is shown in figure [5l where we plot the distribution of triggers in 
(t, p) corresponding to a sine-Gaussian glitch affecting a simplified template bank with 
uniform distribution in total mass. The dependency of SNR and time delay on the 
template mass determines the shape of the cluster. Low-mass templates produce the 
triggers with the largest delay and smallest SNR. The last triggers can be delayed by 
several minutes, which is much longer than the duration of the glitch (typically seconds 
or fractions of a second). 
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Figure 5. Cluster of triggers generated by a template bank responding to a sine- 
Gaussian glitch with to = 0, /o = 20 Hz and Q = 20. The bank consists of 50 3.5PN, 
equal- mass waveforms with M uniformly spaced in [2Mq , 2OM0] . The cluster extends 
well after the duration of the glitch, with the lowest mass templates still producing 
triggers minutes after the glitch. 



5. Conclusion 

Low-mass waveforms represent a potential problem for future inspiral searches and 
in particular for low-latency BNS pipelines such as MBTA [T^ and LLOID [15]. In 
fact, as evident from this study, low-mass triggers generated by glitches could be 
delayed by several minutes. The result is that a single glitch could produce a cluster 
of many spurious triggers lasting for several minutes. Although we found that the 
trigger SNR effectively decreases with increasing delay, glitches can be very strong 
and still produce triggers with a large delay and SNR above the detection threshold. 
Consistency tests such as the x^ test [16] and its variations [IT] are most effective for long 
waveforms and would likely rule out a large fraction of such spurious triggers. Indeed, 
a detailed investigation of the response of these tests is a possible followup of this work. 
Unfortunately however, the efficiency of such tests also decreases with decreasing SNR. 
Triggers with large delay and SNR just above the detection threshold will therefore still 
be problematic and veto procedures based on auxiliary channels will be important. 

We presented three approximations which allow one to predict the SNR and time 
of triggers generated by an inspiral matched ffiter responding to sine-Gaussian glitches. 
We compared them to numerical simulations and investigated their validity in the 
region of (/o, Q, M) parameter space relevant for advanced detectors. Together they 
complement each other, providing full coverage of the explored region. They could be 
useful for determining whether an inspiral trigger can be associated with an excitation of 
auxiliary channels. In fact, the knowledge of the trigger SNR and mass could be used to 
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determine a possible set of (to, /o, Q) parameters and test the hypothesis that the trigger 
was produced by a sine-Gaussian ghtch with such parameters in an auxihary channel. 
Note that sine-Gaussians are the signal basis used by the Omega pipeline [H] — a tool 
widely used by the LIGO collaboration for data quality assessment — and therefore the 
integration of our results into this tool is straightforward. Another possible application 
is searching for clusters of triggers across the template bank whose distribution in SNR 
and time matches the one predicted by the approximations. This could provide a 
way of simultaneously vetoing all triggers belonging to such clusters. The definition 
of procedures for actively exploiting our approximations is another natural followup of 
this paper. 
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